Quantum discord and classical correlations in the bond-charge Hubbard model: quantum phase 
transitions, ODLRO and violation of the monogamy property for discord 
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We study the quantum discord (QD) and the classical correlations (CC) in a reference model for strongly 
correlated electrons, the one-dimensional bond-charge extended Hubbard model. We show that the comparison 
of QD and CC and of their derivatives in the direct and reciprocal lattice allows one to efficiently inspect the 
structure of two-points driven quantum phase transitions (QPTs), discriminating those at which off diagonal 
long range order (ODLRO) is involved. Moreover we observe that QD between pair of sites is a monotonic 
function of ODLRO, thus establishing a direct relation between the latter and two point quantum correlations 
different from the entanglement. The study of the ground state properties allows to show that for a whole class of 
permutation invariant (r;-pair) states quantum discord can violate the monogamy property, both in presence and 
in absence of bipartite entanglement. In the thermodynamic limit, due to the presence of ODLRO, the violation 
for r;-pair states is maximal, while for the purely fermionic ground state is finite. From a general perspective, 
all our results validate the importance of the concepts of QD and CC for the study of critical condensed-matter 
systems. 

PACS numbers: 03.65.Ud, 03.67.-a. 71.10.Fd 
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I. INTRODUCTION 

A very fertile interplay between the theories of quantum- 
information' and condensed matter^ has developed during the 
last decade. On one side, condensed-matter theory has sug- 
gested a wide range of possibilities for the implementation of 
quantum communicationaP and computational'* tasks. On the 
other side, quantum information theory has provided novel 
and deep insights into the physics of condensed-matter sys- 
tems. In particular, since the capability of entanglement to 
mark quantum phase transitions (QPT) was first recognized^, 
the concept of quantum coiTelations has become essential in 
characterizing quantum critical phenomena^ In this con- 
text, quantum correlations and entanglement have been usu- 
ally identified as one and the same concept. It is now accepted 
that the notion of entanglement''' is unfit to capture the whole 
amount of quantum coiTelations present in a system. The in- 
troduction of quantum discord (QD)'"'''-2. in fact showed that 
a system can exhibit quantum correlations even in absence 
of entanglement. Furthermore, QD allows to properly distin- 
guish the total correlations between two subsystems in terms 
of quantum and classicaU^. QD was first devised in the realm 
of quantum-measurement theory'^ but it has been analyzed in 
a large variety of different physical systems. In a quantum in- 
formation perspective, it has been suggested that QD, rather 
than entanglement, may be the most fundamental resource al- 
lowing for the speedups of quantum over classical computa- 
tion'^, and that it may have a relevant role in quantum com- 
munications protocols'^. The study of randomly generated 
states has established that QD is present in almost all quan- 
tum states-". The notion of QD has also been investigated, 
both from the theoretical^'-^- and the experimental side2^, in 
the context of open quantum systems where it was shown that 
QD is generally more robust than entanglement to dissipation 
and thermal noise; moreover, QD is tightly related to complete 
positivity of quantum maps^"*. In the realm of many-body sys- 
tems, the behavior of QD has been analyzed in relation to 



QPTs and thermal effects2^— . So far, the research has mostly 
concentrated on one-dimensional spin 1/2 models 2^2^. The 
main results of these analyses show that two-point QD and 
CC between near as well as distant sites show clear signatures 
of QPTs (discontinuities or divergences), which can be under- 
stood within a general framework^^ and agree with finite-size 
scaling theory in the case of finite chains. Apart from spin sys- 
tems, many-body QD has been studied in the LMG model— 
and the Castelnovo-Chamon modeP', where it marks a topo- 
logical QPT. 

At present, a thorough analysis of QD and CC in correlated 
electron systems is still lacking, mainly because the latter, at 
variance with the simplest spin systems, requires the evalua- 
tion of the discord for pairs of q-dits. The present work is in- 
tended to fill this gap, by investigating the behavior of QD and 
CC for the ground states of the one-dimensional bond-charge 
extended Hubbard model^^'^'', which is a reference model in 
coiTelated-electron theory. The model has an integrable point, 
and its entanglement properties have been the subject of recent 
studies ^'^'^^ where use of two-point and multipartite entangle- 
ment measures led to a classification of QPTs into multipartite 
or two-point driven. These studies left open the problem of ad- 
dressing the general role of bipartite correlations for all two- 
points driven QPTs, as well as their relation with the presence 
of off diagonal long range order (ODLRO) which character- 
izes some ordered phases of the model. The introduction of 
QD and CC allows to solve this problem in a proper frame- 
work. 

We systematically consider the quantum discord and the clas- 
sical correlations, in direct space between two-sites and in mo- 
mentum space between two couples of modes, and we study 
their interplay and their ability to properly describe the rich 
zero temperature phase diagram and the various phase transi- 
tions exhibited by the extended Hubbard model we consider. 
We see that QD and CC can highlight the presence of a so- 
called entanglement transition, where a different role of quan- 
tum and classical correlation at a transtion is revealed by the 
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different behaviour of QD and CC, both in their maxima and 
in the divergence of their derivatives (subsection |IV At . The 
study of the derivatives of QD and CC close to the critical 
lines allows to confirm the two-point/multipartite nature of the 
various transitions and to distinguish transitions that are phys- 
ically different based on a different role of long range quan- 
tum correlations(ODLRO) (subsection llV Cl l. We demonstrate 
that these long-range correlations correlations, which are at 
the basis of superconductivity, are related to two-point discord 
rather than two-point entanglement: indeed a direct relation 
between ODLRO and QD can be found (subsection |IV Bt . 
This relation is true both in the direct and in the reciprocal 
lattice picture, since a functional relation between the two- 
site QD and the two mode QD can be established (subsec- 
tion |IVD]|. As an example of how condensed matter systems 
constitute a natural playground to test quantum information 
concepts, our study of the ground state properties also sheds 
light onto an aspect of the quantum correlations that is very 
relevant in the general context of quantum information the- 
ory: the monogamy property (we address this issue in de- 
tail in the self-contained subsection llV El ). Upon considering 
ground states of the model also at finite system size, we can 
extend previous analyses of the monogamy relation to an n- 
partite setting with n > 3. In a phase of the model the ground 
states coincide with a class of permutation-invariant states, for 
which we show that the monogamy relation is always violated, 
both in presence and in absence of entanglement. In the TDL 
the entanglement vanishes and the violation of the monogamy 
property for QD becomes maximal: due to the presence of 
ODLRO, a single qubit can exhibit finite amount of discord 
with an infinite number of other qubits. The monogamy re- 
lation can be violated also in absence of ODLRO, but in this 
case the violation is not maximal. 

The paper is organised as follows. Section is an intro- 
duction to QD. Section Hn] is a brief review of the main fea- 
tures of the extended Hubbard model: Hamiltonian, phase di- 
agram., etc. Section HV] is the core of the work, displaying 
our analysis of the behavior of QD and CC in the whole phase 
diagram of our model with a special focus on quantum criti- 
cal points/lines. Section [V] closes the paper, highlighting and 
summarizing the main conclusions of our work. 



II. QUANTUM DISCORD AND CLASSICAL 
CORRELATIONS 

The peculiarity of quantum physics mirrors into the highly 
nonclassical nature of correlations between quantum systems. 
It is widely known that quantum entanglement'^ lays at the 
heart of many quantum phenomena and plays a crucial role 
in quantum information processing. However, entanglement 
alone is insufficient to describe the quantum character of the 
correlations present in quantum states. The definition of a 
proper measure of quantum correlations can be derived as the 
difference between the total correlations between two subsys- 
tems A and B, represented by the quantum mutual informa- 
tion I{A : B), and the classical correlations C{A : B), whose 



measure has been introduced ioi^: 

Q{A : B) = I{A : B) - C{A : B) (1) 

where I{A : B) = S{A) + S{B) - S{A : B) is expressed in 
terms of the von Neumann entropy of reduced density matri- 
ces Pa, Pb and the density matrix of the system p respectively. 
The measure for classical correlations is defined as the maxi- 
mal amount of information on one of the subsystems, say A, 
that one can extract by classical means i.e., by operating a 
complete measurement process the other subsystem B 

C{A:B)^m'^x[S{A)-y2pkS{pAk)]- (2) 

k 

Here the set of positive operators {tA®Bk}, 'l2k ^ 
represent a von Neumann measurement, i.e., a set of orthog- 
onal projectors on subsystem B and pk — Ti'Bij^BkpBk] is 
the k-th post-measurement state of subsystem A obtained with 
probability = Ti'APk (see the Appendix for details). The 
maximum over all von Neumann measurements is attained in 
correspondence of the minimum of the conditional entropy 
S{A\B) = J2k PkS{pAk) and in general requires difficult op- 
timization procedures. The definition of Q{A : B) above co- 
incides with the definition of quantum discord originally given 
in—. The reasons that led to this definition and the properties 
of this measure explain how it entails for the quantification of 
quantum correlations. The condition that captures the nature 
of the quantum correlations described by the discord is that if 
a state has non zero discord all complete measurements on a 
subsystem B will unavoidably disturb the whole system and 
in particular the subsystem A: 

p^^lA®BkplA®Bk (3) 

k 

This is a quantum feature that has no classical counterpart, 
but it does not require the presence of entanglement: sepa- 
rable states can have finite amount of discord. The discord 
Q{A : B) — iff the state p is diagonal in a product eigen- 
basis"'^, and thus a classical state has the form {\i)A <8) |j) b} 
i.e., p = K,j\T')AA{i\ «) \j)BB{j\, for a possible classi- 
fication of quantum states see"*". Quantum discord is a corre- 
lation measure aimed at quantifying all quantum correlations 
including entanglement. A key feature of QD is that in gen- 
eral it is not symmetric under exchange of the two subsys- 
tems^: measurement processes applied to the two different 
subsystems can lead to different values of discord. However, 
in the examples we will study the density matrices are sym- 
metric with respect to the exchange of subsystem and so is the 
discord evaluated for the systems represented by these states. 
As already mentioned, the evaluation of QD and CC poses 
in general some difficulties since it requires an optimization 
procedure. For two qubit systems analytical formulas have 
been found only for special examples of two qubit density 
matrices^^, while for continuous variable systems a general 
formula have been derived for Gaussian states only^. In the 
following we will use the analytical formulas^ for the two 
qubit case, while we will resort to a simple numerical opti- 
mization for the case of two qutrits. 
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III. THE BOND-CHARGE EXTENDED HUBBARD 
MODEL 

A. Basics of the model 

The bond-charge extended Hubbard model was derived 
as an effective one-band Hamiltonian for the description of 
cuprate superconductors^^. The model is described by the fol- 
lowing Hamiltonian: 



<i.j>(T ia 

+uY,(n^t-l) (n^i-l) (4) 



where cj^. and Cj^. are fermionic creation and annihilation op- 
erators on a one-dimensional chain of length L; a =t, \. is the 
spin label, a denotes its opposite, n^-^ — Cjcr'^ja- '^he spin-cr 
electron charge, and (i, j) stands for neighboring sites on the 
chain; u and .t (0 < a; < 1) are the (dimensionless) on-site 
Coulomb repulsion and bond-charge interaction parameters; /i 
is the chemical potential, and the corresponding term allows 
for arbitrary filling. 
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FIG. 1: Ground-state phase diagram of H . LEFT: n-u plane; empty 
circles stand for empty sites, slashed and full circles stand for singly 
and doubly occupied sites, respectively. RIGHT: [i,-u plane. 



The model is considered here at x = 1, in which case the 
system becomes integrable. This happens for two reasons. 
First, the 1 — xiriis + njs) term suppresses several hopping 
possibilities. As a result, we can separate the four possible 
states at site i into two groups, namely A ~ { | t) , | i) } and 
^ ~ I ti)}- hopping permutes states of group A with 
states of group B, but not states of the same group. The role of 
spin orientation becomes dynamically irrelevant, and the sys- 
tem behaves as if at each site the local space had dimension 3: 
I t) and I I) can be considered as the same state. Second, the 
hopping term commutes with the terms in u, /i and the num- 
ber of doubly occupied sites becomes therefore a conserved 
quantity. 

The physics of the system described by H is basically that of 
Ns spinless fermions - singly occupied sites - which move in 



a background of L — Ng bosons, of which Nd are doubly oc- 
cupied sites and the remaining are empty sites . Both Ns and 
Nd are conserved quantities, and determine the total number 
of electrons N = Ns + 2Nd. 

The situation may be understood in the formalism developed 
by Sutherland in Ref. — . We can say that, apart from constant 
terms, H acts as a permutator of just two Sutherland species 
(SSs), the Ng fermions, and the L ~ Ng bosons. In practice, 
empty and doubly occupied sites — though different as phys- 
ical species — belong to the same SS, since the off-diagonal 
part of the Hamiltonian (i.e., the hopping term) cannot dis- 
tinguish between them. It is only the constant term counting 
doubly occupied sites which depends on the actual value of 
Nd. 

It is convenient to rewrite both the Hamiltonian and the local 
vector space in terms of the Hubbard-like projection operators 
X°'^ ^ \a),{l3\„ with local algebra Xf^X/* = Sp^X"^ and 
nonlocal (anti-)commutation relations given by 



(-) 



(5) 



here a — 0,1,2, |0)i = |vac)i is the local vacuum, = 
Xl^\0)i is the singly occupied state (with odd parity), and 
\2)i = Xf^\0)i is the doubly occupied state. More precisely, 
as far as the ground state is concerned, the model Hamiltonian 
in the one-dimensional case can be fruitfully written as 

H = - Y.{Xl'^X^l,-Xf'Xll,+H.c.)+uY^Xf^ 

i i 

' (M+^)E(^"+2^n ■ (6) 

i 

The eigenstates are easily worked out^^*^, and read 

Nd) AA(ryt)A'.lio . . . ^lo^ |vac) ; (7) 

Here Af = [{L - Ns - Nd)\/{L - NsY-Nd^]"'"^ is a normal- 
ization factor; X^^ is the Fourier transform of the Hubbard 
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projection operator i.e., X^^ = -^exp{ij^jk)Xj . 

Moreover, rj^ = -^i'^ ^1^° known as the eta operator, 

commuting with H; (77^^)^'* creates Nd pairs which are fully 
spread over the chain. These are the 77 pairs first introduced by 
Yang"*^. This structure corresponds to a very simple physical 
picture: eigenstates contain Ns spinless fermions in momen- 
tum eigenstates {^^^, . . . , HZI^^X^j ^nd Nd spinless bosons (rj 
pairs). 

The energy eigenvalues are given by E{{nk},Nd) = 
2fiNd - (m + ^)Ns where rik = 0, 1 



is the number of fermions with momentum For any 

given Ns = J2k and Nd the minimum is achieved by oc- 
cupying with Ns fermionic particles the momentum modes 
{7r(A^s — 1)/L, . . . , tt{Ns — 1) /L}, the corresponding eigen- 
value being £:(7Vs, A^d) = -2 sin (tt-^) / sin {j^) ~ 2 ^iNd- 
(fi + Ns , whence we obtain the ground state energy den- 
sity in the TDL 



2 / 
£{ns,nd) = sin(7rns) - 2fmd - ( M + 

TT \ 



(8) 
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with £ = E/L, Us = Ns/L, Ud = Nd/L. 

The actual ground state is found by requiring that Ug and rid 

minimize ([8]). 

For /J, < we have Ud = 0, hence upon minimizing we get 
Us = ^ arccos (-f - |). For -4 - 2// < u < 4 - 2/i 
we have empty and singly occupied sites (phase I), for u > 
4 — 2/x we have only singly occupied sites (phase IV) and for 
M<— 4 — 2/i we have only empty sites. 
For /X > we have Ud — (1 — ris), hence upon minimizing 
we get ris ^ arccos - |) . For -4 + 2/i < u < 4 + 2/i 
we have doubly and singly occupied sites (phase I'), for 
u > 4 — 2/i we have only singly occupied sites (phase IV) 
and for u < —4 — 2/i we have only doubly occupied sites. 
For /J, = we get Us = ^ arccos (— f )• For —4 < u < 4 we 
have empty, doubly and singly occupied sites (phase II), while 
for u > 4 — 2/t we have only singly occupied sites (phase IV) 
and for u < —4 — 2/i we have empty and doubly occupied 
sites (phase III). 

Hence we get in the /i-u plane the phase diagram depicted in 
Fig. [T] right. In the left part, the same ground-state phase 
diagram is drawn in the n-u plane (with n = N/L average 
per-site filling). The phase diagram presents various QPTs 
driven by parameters u and /t (or n). Each transition is char- 
acterized by a change in the number of on-site levels involved 
in the state. Phase IV has just one level per site since each site 
is singly occupied. Phases I and I' (which is the particle-hole 
counterpart of phase I) have two on-site levels: singly occu- 
pied sites and empty or doubly occupied sites respectively. 
This holds for phase III as well, where only empty and doubly 
occupied sites appear Phase II is the only phase in which all 
three on-site levels are involved. Phases II and III are char- 
acterized by the occurrence of off-diagonal long-range order 
(ODLRO) and superconducting correlations, evaluated as; 

lim {Xf°X^^^) = udil -Ud- Us) . (9) 

r— f oo 

Note that ODLRO — though not allowing real superconduct- 
ing order at x = 1 due to spin degeneracy, which implies the 
vanishing of spin gap, is at the very root of superconducting 
order, which occurs atxy^ 1.*'^ 

Before discussing the various transitions in terms of the QD 
behavior, let us recall some feature of each of them in terms 
of standard theory. First of all, since Nd and Ns are both con- 
served quantities, the transitions should be originated from 
level crossing. Indeed, they also occur at finite L. Never- 
theless, none of them is of first order, since it can be easily 
checked that the first derivative of Egs is always smooth. In 
fact, the transitions I (F) ^- IV and II IV and II III 
are second-order QPTs, while the transition II I (F) is an 
infinite-order QPT. 



B. Reduced density matrices 

The present work focuses on two-point correlations. To 
evaluate them, knowledge of the ground-state reduced den- 
sity matrices is necessary, and we shall report their expression 
for completess (for a full derivation, the reader may refer to 



Refs. • ). Correlations can be analyzed within two differ- 
ent and complementary pictures. Obviously, one can exam- 
ine correlations between sites of the lattice (direct lattice pic- 
ture). In addition, the structure of eigenstates in the model 
suggests yet another approach, namely to consider the recip- 
rocal lattice, whose elementary nodes are momentum modes 
kj = ^j, j — 0, . . . , L — 1. In some respects, the reciprocal 
lattice picture affords a simpler description of the system.^'*'''^. 
Let us start by giving reduced density matrices the direct lat- 
tice picture. The one-site reduced density matrix pi when ex- 
pressed in terms of the basis {|0), |2)},i is diagonal in all 
the regions of the phase diagram: 

Pi ^ diag {l-Hs-Ud, Ud} , (10) 

the two-site reduced density matrix pij in the basis 
{|00), |01), |02), |10), 111), |12), |20), |21), |22)},, reads^: 
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with c = n<j/(l - Us), Pij = (1 - Tis)^ - |rijP. and 

-p _ sin(ws7r|i-j|) 
7r|«-i|) • 

Let us now turn to the reciprocal lattice picture. To each mo- 
mentum mode kj corresponds a 4-dimensional Hilbert space, 
spanned by the basis 

Bfc, = io>fc,, it>fc,, u>fc,, \n)k„ (12) 

The reduced density matrix for any such mode reads, in the 
TDL, and in the basis (fT2l i. 

Pk^ = diag{a'^,ab,ab,b'^) (13) 
where a ^ and b = -r^. 

The two-mode (16 x 16) reduced density matrix for modes ki 
and kj, ki ^ kj, is diagonal with respect to the local basis 
®Bk -- In the TDL, the eigenvalues are a"b*^" with mul- 
tiphcity TTia = (l) ■ 

The case ki ~ ~kj has to be treated separately. The two- 
mode (16 X 16) reduced density matrix for modes kj and —kj 
has support on a 4 x 4 subblock. Indeed the sole states that 
can be built by the action of the 77^ operators belong to the 
subspace spanned by: 

6fc^,_fe^ = {|o,o)j, I t,;>j, I ;,t>j, I n,n>j} (i4) 
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where \a f3)j = \a)kj «) |/3)-fe,-- 

In the TDL, and in this basis, the nonvanishing subblock of 
the matrix reads: 



Pie, 





ab ab 

ab ab 

V fe2 



C. Behavior of entanglement at QPTs 

Two-point entanglement at the QPTs of the model was thor- 
oughly analysed in Refsi^^, upon consideration of different 
correlation measures: the two-point concurrence Ki j or the 
two-point negativity Afij as measures of entanglement (notice 
that definition of concurrence is available for two-qutrit sys- 
tems), the mutual information li j as a measure of total two- 
point correlations, and the single site entropy Si as a measure 
of multipartite entanglement between one site and the rest of 
the chain. The behavior of all correlation measures was stud- 
ied as a function of x (x — fi or x ~ u) in the vicinity of the 
quantum critical points. Results are briefly summarized in the 
table below. 



out in the Appendix, see lA 11 1. Evaluating the mim- 
imum of the reduced conditional entropy reduces, Eq. 
(IA7l l. to taking the minimum among two functions, i.e., 
in^Bk}S{pij\{Bk}) = min{5'i,5'2}, where 5*1, 5*2 depend 

on 01 = ^(l + 4n^-4n)+4|r,,f ,02= , 
(15) ^ |i+2„^-3n-2|r..n g3j|A6j. 

We verify that for all values of |i — j| we always have 5*1 < S2 
and therefore two-point classical correlation ^ and quantum 
discord ([T]i can be written in terms of . 
In order to compare quantum discord and entanglement, we 
also evaluate two-point concurrence^^ 
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finite 



finite 



finite 
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Q2 



QS 



Q2 
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The analysis of divergences allows to classify the different 
transitions into those driven by two-point correlations (Q2: 
II — > III, II — IV, I — >^ II), where some two-point correla- 
tion measure {Ki,j, Afij or diverges, and those driven 
by multipartite correlations (QS: II — > I) where only Si di- 
verges. However, the two-point charachter of the transitions 
II III, II — T' IV is only detected by /, j (a measure of total 
correlations), while Ni,j (the measure of quantum correlations 
used) is unfit to discriminate between those transitions and the 
multipartite-driven one (II — > I). 



IV. RESULTS 
A. Region I (I') 



K,,, = mill {0, |r,, - \/((i-^)'-ir,,f)(n2-|r.,f)| } . 

(16) 

In the following, letters Qij,Cij,Ii,j,Kij alway denote 
quantum discord, classical correlations, mutual information 
and concurrence respectively. The values of lij, Cij, Qij 
and Kij for region I and different values of | z — j | are plotted 
in Fig. 121 
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FIG. 2: quantum mutual information (blue, solid), quantum dis- 
cord Qi,j (red, solid), classical correlation dj (green, solid), con- 
currence Kij (red, dashed) as a function of fi in region I (u = 4), 
for — ~ 1 (top, left), — = 2 (top, right), — ~ 3 (bottom, 
left), |j - j| = 4 (bottom, right). 



We start our analysis by evaluating correlations (QD and 
CC) in phase I. Results for phase I' are omitted, since they are 
exactly equal (by virtue of the particle-hole symmetry one just 
has to replace emply with doubly occupied sites). 
Phase I (T) is characterized by the absence of doubly occupied 
(empty) sites, so that the effective number of on-site levels re- 
duces to 2. Consequently, the 2-site reduced 9x9 density ma- 
trix pij has nonzero entries only in the 4x4 subblock spanned 
by{|00),|01>,|10>,|ll>}.,. 

The quantum discord can be evaluated analytically through 
the methods developed in Refs"^^^^ (details are worked 



We first see that the Qij and C'ij have the typical oscillat- 
ing behavior already shown by the mutual information^^. At 
variance with the previous analysis, where the quantum cor- 
relations measured by the concurrence were different from 
zero only in proximity of the borders of the regions i.e., for 
/i — > —4, 0, here we see that the system exhibits non zero dis- 
cord within the whole region I except at some nodal points de- 
fined by the equation — {cjcj) — '^'"^Jj°^|^'|y-'^'' = where 
all correlation measures vanish, J — Cij = Qij = 0. Clas- 
sical correlations show a similar behavior Therefore, in the 
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central region of phase I, where Kij vanishes V|j — j\ > 1, 
two-point discord and classical correlations are still present. 
Correlations are modulated by the sinusoidal behavior in- 
duced by Tij and at fixed /i they all decay algebraically with 
the distance: lij, Qij, Cij ~ |z — see Fig|3] 



I, C. Qatfixe(lp=-0.1 
0.1 




FIG. 3: quantum mutual information (blue, solid), quantum 
discord Qij (red, solid), classical correlation dj (green, solid) as a 
function of |i — j| in region I for /x = —0.1, u = A. Upper dashed 
lines represent the envelope of the respective maxima wich exhibits 
a power law decay (~ | i — j | ~ ^ ) 



In proximity of the transition I — IV it was shown in'''' that 
the system exhibits an entanglement transition^: the range 
of the entanglement TZjc, i.e., the maximal distance \i — j\ 
for which Kij ^ 0, goes to infinity when approaching the 
transition. In particular, Ki^j have a maximum value for 
TT-s — )■ 1 as |i — j| oo. This behavior is reflected in that 
of /i j, Qi.j, Ci,j, which also exhibit a global maximum at a 

value ni*'-''' « 1 — l/(2|i — j\) which approaches Ug = 1 for 
|i — j| — > oo. Hence, the behavior of discord miiTors that of 
the entanglement. This behavior is depicted in Fig|4] In fact, 
also the mutual information and the classical correlations ex- 
hibit the same kind of behavior. However the values of the 
maxima for the various measures , , Qij , Cij scale in a dif- 



ferent way with the distance: 



Q.K 




FIG. 4: maxima of quantum discord Qij (solid lines) and concur- 
rence Ki,j (dashed lines) for \i — j\ = 16 (blue), |i — i| =32 (red), 
|i - j| = 64 (green), 
region l(u — 4) 



- j\ ~ 128 (black) as a function of in 



Iin^:''\\^-J\ 
Cini^'^\\t-j\ 



1 



log(l»- j|) 



(17) 

(18) 

(19) 
(20) 



Therefore, when approaching the metal insulator transition 
I^IV the maxima of correlation measures {lij^Qij^Cij^ 
decay in algebraic way along the chain. Quantum discord 
dominates for high distances, since the spreading of the clas- 
sical correlation is suppressed by a factor log(|« — — j\- 
This difference in the behavior of QD and CC defines the 
different role that they have at the transition and can be fur- 
ther appreciated by studying the derivatives of the different 
coiTelation measures with respect to /i. In the critical limit 
/i — 0, —4 we have 



while 



d^iC{pij) ~ log l/i - pc 



(21) 



(22) 



Therefore, while the d^Qi j correctly agrees with the scaling 
behavior of d^Ii^j and d^Ki j evaluated in-'^, df^Cij though 
being singular has a lower degree of divergence, so that clas- 
sical correlations are subleading in the vicinity of the critical 
point. 

We therefore see that the introduction of the new measures of 
correlations Qi j and Ci,j and the study of their derivatives 
allows on one hand to properly identify the metal-insulator 
transition and to properly classify it as a two-point QPT^», 
and on the other handallows for a refinement in description of 
the QPT. The importance of this feature will be more evident 
in the following paragraphs where we will describe the other 
two-point QPTs i.e., II IV and II— II. We close this subsec- 
tion by discussing the role of the divergences of the different 
correlation measures and their relation with the divergences of 
the energy density of the system. In " , the authors found a di- 
rect relationship between the singularities (discontinuities and 
divergences) in the derivatives of the energy density of the sys- 
tem £ = E/L with respect to the parameter A that drives the 
QPTs, and the singularities in the elements of the two-point 
reduced density matrix pij or their derivatives with respect to 
A. In our case the, the divergences in d\Iij and d\Qi,j inherit 
the non analyticities of the derivatives of the elements of pij 
at the critical point. In particular the elements: 



dxD2,dxOi,dx02 -> 



1 



(23) 



show the same divergences exhibited by the second deriva- 
tive of the energy density dHJ with respect to X = p (I— >IV), 
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FIG. 5: quantum mutual information I (blue), quantum discord Q 
(red), classical correlation C (green) as a function of Ud in region III 



i.e., dfi^S ^1/ ^y\fJ~~JL^. However, as we have seen above, 
classical correlations, though diverging, show a logarithmic 
divergence instead of an algebraic one, and accordingly one 
might believe that this is an accidental fact due to the defini- 
tion of the correlation measure (i.e., that CC always display a 
lower degree of divergence). However, as we will see in the 
next sections, the classical correlations Cij behave like Qij 
and lij, in terms of their derivatives with respect to A = u, at 
the transition II— and therefore they coherently behave as 
the energy density at that transition i.e., d^£ ~ l/y/\u — Uc\. 
In summary, while the derivatives of different elements of pij 
and of some of the correlations measures defined on the pij 
show the same divergent behavior at the various transitions, 
which agrees with that of the energy density the classi- 
cal correlations may show different kind of divergences and 
are thus able to discriminate between quantum phase transi- 
tions that are physically different. 



tangled correlations but rather in the form of QD. Moreover, 
we notice that QD, as well as CC, between any two sites has 
the same value, irrespective of their distance: this reflects the 
way in which the 77-paring mechanism spreads the correla- 
tions equally along the whole chain. The r^-pairing is also the 
ground for the appearance of ODLRO, which follow directly 
from (|9]l. It is intuitive to suppose that these superconducting 
correlations might be related to some kind of two-point quan- 
tum correlations, and indeed many authors have tried to find 
such a relation, see for example^'''^^ While a relation with 
the entanglement properties in k space was found in^*^ in the 
case of for 77-pairs and BCS states, in direct space this relation 
could not be established in terms of the concurrence since the 
latter vanishes in the TDL"*^. 

While ODLRO in 77-pair states cannot be related to two- 
point entanglement, our analysis allows instead to connect the 
ODLRO to the two-point quantum discord. Indeed we find 
that in the TDL Qi,j = Q^"^^, Vi, j and we have 

QTDL ^ y^^) ^ arctanh(l - 2x) + a; log 16 -f-24) 

log(-l - -1+V1-4J + log( (^^) + log^] 

where f{x) a monotonically increasing function of x = 
nd{l — rid), i e-, of the ODLRO. The above analysis allows 
to establish a direct relation between a fundamental quantum 
property such as ODLRO and the presence of two-point ( two- 
qubit) discord. It therefore seems that the important two-point 
quantum correlations necessary in direct space for the appear- 
ance of the ODLRO are represented by the discord and not by 
the entanglement. 

We finally note that the presence of the ODLRO in 77-paris 
states is reflected also by the behavior of CC, which also are a 
monotonically increasing function of rtd(l — n^). The relation 
between CC and ODLRO will be important in the discussion 
of the transitions described in the next section. 



B. Region III: discord and ODLRO 

Phase III is characterized by the absence of singly occu- 
pied sites, so that the number of on-site levels effectively re- 
duces to 2, and the quantum discord can be evaluated analyt- 
ically in the same way as above. Moreover in this case the 
number of Sutherland species reduces to 1 . The quantum dis- 
cord may be evaluated in the same way as above. We have 
mf[Bk}SiPij\{Bk}) = minis'!, S'2}, where 81,82 depend 
on 6*1 = + and 52 on 6*2 = 1- 2nd, 6*3 = 1-2^^ 

(|A3l|A6l l. Since two-site density matrices pij are equal for all 
the values of two-site correlations are equal for ech pair 
of sites, lij = I,Cij = C,Qij = Q. We have 81 < 82 and 
therefore the classical correlations (|2|l and the discord ([T]i can 
be written in terms of 81 . 

The values of /, C and Q for region III are plotted in Fig. |5] 
The first result of our analysis is that while in the TDL the 
concurrence Kij = min{0,— 2n^(l — Ud)^} — vanishes 
everywhere in region III, the discord is always different from 
zero in the region; we thus have that the ?7-pair states display 
two-point quantum correlations, though not in the form of en- 



C. Region II 

Region II contains empty as well as singly and doubly oc- 
cupied sites, so that there are 3 on-site levels. This means that 
the evaluation of discord and classical correlations is more dif- 
ficult than in the previous cases. In order to evaluate Qi j and 
Ci,j we used the two numerical recipes described in IA 21 The 
two methods show perfect agreement in the value of the dis- 
cord throughout the whole region, and this is a first indica- 
tion of their reliability. A further element of confidence in the 
methods used is the fact that Qij and Cij must be continuous 
in the transitions II 1 , III (since all matrix elements of pij 
are): when we approach the phase boundaries, the numerical 
limits of Qi j and Cij in region II coincide with the analytical 
values determined in region I and III. 

In Fig. |6]and Fig. |2]we plot li j, Ctj, Qij in region II 
as a function of u for |« — j| = 1, 2, 3, 4, and for n = 1 and 
n — 0.5 respectively. In the Table below, we summarize the 
critical behavior of the derivatives of quantum discord Qi j 
and classical correlations Ci,j for the transitions II I, II 
— > III, IV. These values are obtained as follows. We find 
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FIG. 6: quantum mutual information (blue), quantum discord 
Qi^j (red), classical correlation Ci.j (green) as a function of u in 
region II for n = — j\ = 1 (top, left), \i — j| =2 (top, right), 
|z — j| =3 (bottom, left), |i — j| =4 (bottom, right) 




-4.0-3.5-3.0-2.5-2.0-1.5-1.0-0.5 -4.0-3.5-3.0-2.5-2.0-1.5-1.0-0.5 



FIG. 7: quantum mutual information (blue), quantum discord 
Qij (red), classical correlation Ci,j (green) as a function of u in 
region II for n = 1, |i — j| = 1 (top, left), \i — j| =2 (top, right), 
j — j| =3 (bottom, left), |i — j| =4 (bottom, right). We note that at 
the multipartite transition II— >I (u = 0) all the two-point correlation 
measures behave in a smooth way 



numerically (with either of the procedures sketched above) 
the optimal measurement which minimizes the reduced 
conditional entropy. Contrary to what happens in region I, the 
orthogonal measurement minimizing the conditional entropy 
varies throughout region II, i.e., the parameters of the unitary 
rotation V are not constant throughout the whole region. 
However, in the neighborhood of the critical lines {u —4 
and u — > — 4cos7rn) they are found to remain constant at 
any fixed n. We therefore use these constant values in the 
expressions for the reduced conditional entropy and obtain 
analytical formulas for Qi^ and as a function of u. 
We then extrapolate the critical behavior by studying these 
functions in the critical limit. 



transition 



11^ I 



II III 



11^ IV 



II III 



1/2 



1/2 



^ 4 



^ -4 



dl_ 

All 



finite 



dQ 



finite 



dC 



finite 



oc log(uc - u) 



The results can be summarized as follows. 



II — > I o (or I') two-point Qij, hj, C, 



In the transition 
are regular, thus 
confirming that this transition has a multipartite nature. 

As for the transitions II^IV and II^III, previous 
analyses^^ have shown that both transitions have a two-point 
character As a first result, we see that at both transitions 
quantum discord is able to correctly detect the divergence 
expected, whereas negativity fails for this aspect^"* (see Sec. 
unci ). The two transitions are however physically inequiv- 
alent, since they lead to two completely different phases: 
transition II-^-IV is characterized by the disappearance of 
ODLRO, whereas at transition II^III ODLRO is present. We 
now show how this difference can be properly described by 
the study of the two-point classical correlations. 
In the transition II -> III, while duIi.j,duQi,j > and 
duCij < all the derivatives display the same kind of alge- 
braic singularity. On the other hand, in the transition II — > IV, 
we have that duhj, duQi.j, duCij < 0, they all diverge, but 
Cij has a lower degree of divergence i.e., it is logarithmic; 
this property allows to correctly describe the transition as a 
two-point one and furthermore to assimilate it to the metal- 
insulator transition I— >IV where the CC show the same kind 
of divergence. 

The result can be further deepened by considering the fol- 
lowing argument. All two-point correlations in region II 
can always be split into a finite and an infinite range con- 
tributions: A(j ^ A.^j + A'J^, where A = C,I,Q and 



The infinite range contributions can be analytically evaluated 
and they all explicitly depend on the value of the ODLRO 
in this phase, nd{l — Ug — rid)- Therefore, their derivatives 
with respect to u have all the same behavior: they display the 
same type of algebraic singularity in case of transition II— ^ III 
(ODLRO), while they do not display any singularity in case 
of transition II— >IV (disappearance of ODLRO). 
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On the other hand, as for the finite range contributions we 
find that d^Cij diverges at both transitions but with a log- 
arithmic behaviour that is dominant only in the transition 
II— >IV (where 9„Coo is regular) while its quantum counter- 
part duQi,j diverges algebraically. The above results show 
that the introduction of the discord and classical correlations 
allows to discriminate between two apparently similar but in- 
equivalent two-point QPTs, and to root their difference in the 
persistence (disappearance) of ODLRO at the transitions. 



D. Reciprocal Lattice 

We now consider quantum discord between two momen- 
tum modes in the reciprocal lattice; the analysis is significant 
in region II and III, where 77-pairs are present, and for val- 
ues of kj > kg where kg = 27^^ maximum single- 
fermion momentum, since the portion of k-space pertaining to 
single fermions is factorized. Let us fist consider two modes 
kj kj. From the results derived in^*" we have that the 
measures of correlations all depend on a single parameter a 
linked to the average occupation number of a generic mode 
kj, a — (nfc^-)/2 = n^/ {1 — ris)^ Vfcj. In particular, the 
only pairs of modes (fc^, kj) which are correlated are the ones 
for which ki — —kj, while if ki ^ —kj the relative mo- 
mentum modes are completely uncorrected i.e., Ik^.kj = 
and therefore Qki.kj — 0. When ki — —kj the single- 
mode von Neumann entropy reads S{pkj) = — 2(aloga + 
(1 — a) log(l — a)), the two-mode von Neumann entropy is 
S{pki,kj ) = S{pkj ) + 2a(l — a) and hence the mutual infor- 
mation is Iki,kj = — 2(a log a + b log b — ab). 
In order to evaluate the quantum discord, we should now con- 
sider the reduced conditional entropy after a generic measure- 
ment is performed on mode kj, and minimize with respect 
to all measurements. It turns out that, a von-Neumann mea- 
surement B = {Ho, n^, n^, n^^} onto the trivial basis B-kj 
yields 

Po = — Tr_fc^.nopfc^,-fe^no = a2|0)(0| 
Po 

Pt = —Tr^k,U^Pk,,-k,U^=ab\t){t\ 
Pi = ^Tr_fe^n;pfc^,_fe^n^ = ab\ I 
PU = —Ti-k,U^iPk,,-k,U^i = b^\n){n\ (25) 

Pn 

so that J2aPaS{pa) = and the minimum is immediately 
attained. Consequently we have that the quantum discord has 
a simple expression 

Qk,-k, = Ik,-k, - S{pk^) = 2a(l - a) cx Nk^-k, (26) 

and it is simply proportional to the negativity Nkj-kf'''- ■ 
This result allows us to derive some important conclusions. 
On one hand, the relationship found in^^ between ODLRO 
and negativity in region II can be rewritten in terms of the 



l.Or R 




FIG. 8: Left: ratio R in phase III for chains of varying length L, with 
Nd = l/L (red), A''^ = [L/2J/L (green). Right: ratio R in phase I 
in the TDL. 

discord Qk .-k showing once again the quantum roots of the 
ODLRO: 

lim (Xf ) = {l-ns)^a{l-a) = {l-n,fQk,,^kj2. 

(27) 

This result, together with eq. |25] allows to establish a func- 
tional relation between the two-site discord Qi j and the two 
mode discord Qk^-ki- 

On the other hand, the line n = 1 is an iso-correlation 
line-'^: since a = 1/2 = const, and therefore the momentum 
particle density {ukj), and all the correlations between sub- 
system in the momentum picture are maximal and constant 
in the whole phase II. In particular, they are constant at the 
transition II— >III, therefore this transition cannot be identified 
by studying the derivatives of the correlation measures in k 
space. On the other hand, at the transition II^-IV there is a 
sudden change in all correlations that discontinuously drop to 
zero in correspondence of the insulating phase that is charac- 
terized by a ground state which is factorized also in the mo- 
mentum space i.e., = ®kj |<7)fcj^ with a =t,i- We there- 
fore see that the two-point transitions II— > III and II— >IV can 
be distinguished even in momentum space, and this reinforces 
the result obtained in the previous section in the direct lattice 
picture, where the difference between the two transitions is 
highlighted by the behavior of Ci^j. 

E. Monogamy of quantum discord 

The study of the ground state properties of the extended 
Hubbard model can be fruitfully extended in order to assess a 
relevant quantum information problem: the relations between 
entanglement and discord. In this framework one of the main 
question to address is the possibility that the discord may sat- 
isfy a monogamy relation similar to the one satisfied by the 
squared concurrence in the case of n-partite qubit states^^: 

^l,(2...L)>^12 + ---^?L (28) 

where if i_(2...l) is the concurrence between site 1 and the rest 
of the system while Ki,j is the concurrence between qubit 1 
and j. Indeed, an analogous relation for the two-qubit discord 
has been recently discussed for the limited case of tripartite 
states case in^^ where it was found that the discord can be 
polygamous for tripartite W states. Here the discussion can 
be extended for the case of more general n-partite pure states. 
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with n > 3; in the following we will focus on two different 
case: the ground states in region 111 and 1. 

As for Region 111, we have that the ry-pairs states coincide 
with class of two-qubits permutational invariant states that can 
be written as 

|^(iVd, L))=( ^] J2 P\Nd, L - Nd) (29) 



where is the binomial coefficient, IN^.L — N^) is 

a fixed state with a given sequence of Nd ones (pairs) and 
L — Nd zeros (empty sites), and the sum is taken over all pos- 
sible permutations Ps (the 3-partite W state belongs to this 
class of states). For states of these form, at fixed Nd < L 
the single and two site reduced density matrices can be easily 
evaluated from ( fTOl l and ( fTTT l and they are equal for all sites, 
i.e.. Pi = pi and pij ~ p2 and the discord can be evaluated as 
described in the previous sections. Since \tp{Nd, L)) is a pure 
state the QD between one site and the rest of the chain is equal 
for all sites Q{i\2.---.l) — Qi ^"d it simply coincides with the 
entanglement between the site and the rest of the chain, i.e., 
Qi = S{pi), it is a function of Ud only and and it is bounded 
by 1. Both for finite Nd,L and in the TDL the two point 
QD Qi j — Q-ii^diL) does not depend on j and therefore 
Yjj Q2{nd, L) = {L ~ l)Q2("-d, L)- As already mentioned, 
similar arguments can be applied to the concurrence K: with 
TT-s = the dependence on \i — i\ disappears and, in particular 
for large L one has Kij w 1/L; for finite Nd, L the con- 
currence is small but different from zero, and the monogamy 
property is always satisfied by the squared concurrence. 
On the other hand, a direct evaluation of the above quanti- 
ties shows that VNd and L > 3 one has R — Qi/[{L — 
^)Q2{nd, L)] < 1; in Fig. [8] (left panel) we show the ratio 
R for Nd = l /L, [L/2\ /L and different values of L. 
While a general analytical demonstration of this result is not 
straightforward, one can note that in the case of permuta- 
tional invariant states, for any fixed value of nd it is always 
possible to find an infinite number of states \'4>{Nd, L)) with 
L = Nd/rid and such that Qi < (L — l)Q2{ndTL), i.e., 
the monogamy relation is violated. Indeed, while Qi just 
depends on nd, Q2{Nd,L) is a decreasing function of L 
which is lower bounded by its TDL expression dZST l. There- 
fore all the states for which L — Nd/nd satisfies the relation 
Qi ^ {L — 1)Q2^'" will violate the monogamy relation. 
As for the TDL, while Ki j — > 0, Qg is constant at fixed nd 
and Q2 = as in and therefore R ^ 0. 

Since the above arguments apply to a whole class of permu- 
tation invariant n-partite two-qubit states (|29] l, we can state 
in full generality a property of two-qubit QD: for n-partite 
states (n > 3) QD can be polygamous both in presence (for 
finite Nd, L) and in absence (TDL) of two point entanglement. 

While it is tempting to relate the violation of a monogamy 
relation by the discord to the presence of those correlations 
that are typical of 7]-pairs states, and that give rise in the 
TDL to ODLRO, our model shows that there are other classes 
of states in which such violation can occur. Indeed, in 
Fig. [Hl^right panel) we report the ratio R for the ground state 



of region 1, which reads: 

mNs,L)) = |fci, . . . , fcjvj = • • • Xll^ I vac) (30) 

i.e., contains Ns fermions in momentum eigen- 
modes (ki,...,kNj, created by action of the 
Fourier transform of the Hubbard projection operator 
^k" = ^exp(ifjfc)Xjo onto the vacuum. The 
results refer to the TDL case and they show that for such 
states, although Qi j does depend on the distance \i — j\, 
the monogamy property is violated by the two point discord 
in proximity of the QPT 1 IV. This feature reflects the 
fact discussed in Section IIV Al that when /i — > 0, there is a 
spreading of the quantum correlations over the whole chain. 
Indeed, the violation of the monogamy condition starts in 
correspondence of w —0.2, where the entanglement has 
already started to spread along the chain and has a finite range 
(TZn diverges only at the transition). 

This result has two interesting consequences. On one hand 
the ground states in region 1 show that, depending on the pa- 
rameters that define them (ng in this case) for the same class 
of states the discord may or may not violate a monogamy 
relation^^^. On the other hand the behavior of the discord 
allows to refine the description about region 1 carried out in^-. 
There, by means of the entanglement and correlation ratio 
it was pointed out that the ground states in region 1 have a 
truly multipartite character in the center of the region, while 
when approaching the transition the weight of the two-point 
correlations starts to increase; and this agrees with the two- 
point character of this transition. Here this picture is revealed 
by the violation of the monogamy property displayed by the 
QD: in order to prepare the two-point transition at p — 0, 
the system reorganizes its correlations in such a way that 
their two point character starts to prevail; one can therefore 
identify the point in which this process starts with the value 
of the parameters i.e., p w —0.2 at which the monogamy 
property is violated by the discord. 

We finally compare the two above cases in terms of the vio- 
lation of the monogamy property. Here the key observation is 
the different kind of violation exhibited by the discord. In re- 
gion 1 the discord can be polygamous but the amount of quan- 
tum correlations shared by a single site with the other sites of 
the chain is finite i.e., < i? < 1 for p ^ and it vanishes at 
the transition p = because Qi 0, while J^j Qi-j tends to 
a finite value . On the contrary, for 77-pair states the violation 
has a completely different nature: each site can be equally 
correlated with all the other sites of the chain: R = OVrtrf. 
This difference is indeed rooted in the presence of ODLRO in 
the TDL and in the previously found relation between discord 
and ODLRO. This kind of violation is associated by the dis- 
appearance of the two-sites entanglement, while for the state 
in region 1, the violation occurs in presence of bipartite entan- 
glement. 

The above results allows to give a general statement about 
quantum discord for multipartite pure states: it can be non- 
monogamous both in presence and in absence of bipartite en- 
tanglement. However the violation of the monogamy property 
can be maximal when ODRLO is established in the TDL and 
no bipartite entanglement is present in the state. 
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V. CONCLUSIONS 

In this paper we have addressed several important questions 
related to the ground state correlation properties of a refer- 
ence fermionic model, the bond-charge Hubbard model. We 
have applied the recently developed measure of quantum dis- 
cord (QD) and classical correlations (CC) to study how these 
relate to quantum phase transitions displayed by the model. 
By means of analytical and numerical analysis we have de- 
rived and analyzed the expressions of QD and CC for two- 
qubits and two-qutrits systems both in the direct lattice and 
in momentum space. Our results allow to describe the differ- 
ent quantum phase transitions in terms of the divergences of 
the various correlation measures. As shown in^^"^^ the tran- 
sitions can be classified on the basis of the relevance of the 
two-point and multipartite correlations involved. At variance 
with other entanglement measures^'*, such as negativity, QD 
(and CC) exhibits the expected non analyticities that define 
the two-point transitions. Moreover, the comparison of their 
behaviour allows to discriminate between two apparently sim- 
ilar kind of two-point transitions. In particular, a careful study 
of the contributions in which CC can be decomposed gives 
the possibility to detect the presence (disappearance) of the 
off diagonal long range order (ODLRO) and to identify its 
consequences at the various transitions. 
Furthermore, the study of the discord between two generic 
sites i, j and two momenta modes ki,—ki allows to establish 
a direct relation between ODLRO and the two-site/momenta 
modes discord, which turns out to be a monotonic function 
of ODLRO. This result is remarkable, since in the TDL no 
two-site entanglement is present in this states. By means of 
the same analysis it is possible to establish a functional re- 
lation between the two-sites discord in direct space Qi j and 
two-modes discord in momentum space Qki-ki- 

The study of Qi j for 77-pairs states is also important for 
describing the behavior of the discord with respect to the 
monogamy propertyiS. Indeed, the 77-pairs states are isomor- 
phic to a relevant class of permutational invariant multipartite 
qubit states. While in the finite size case, all the states in the 
class display non-zero two-qubit entanglement, in the TDL 
the latter vanishes. However, in both cases we have shown 
that two-qubit discord is in general different from zero and 
furthermore it violates a monogamy relation. Finally, we have 
shown for another class of states, the non-interacting fermionc 
ground states in region 1, the discord can be polygamous de- 
pending on the values of the parameters. The main difference 
between the two class of states analyzed resides in the kind 
of violation of the monogamy property: only for the yy-pair 
states the single qubit can be arbitrarily correlated with all 
the other infinite sites, thus leading to a maximal violation 
of the monogamy property. This fact is rooted in the presence 
of ODLRO in these states and in the direct relation between 
ODLRO and discord. 

Our results confirm that the application of quantum infor- 
mation concepts to condensed matter systems can fruitfully 
lead to a precise description of the role of correlations in quan- 
tum phase transitions and at the same time to the development 
of useful relations that shed new light on the nature of quan- 



tum correlations as measured by discord. 

Appendix A: Evaluation of discord 

The difficult step in evaluating the discord Q is the mini- 
mization of the conditional entropy S{pij\{Bk]) with respect 
the set of all von Neumann measurements. 
In general, the concavity of the conditional entropy implies 
that its minimum is attained through extremal POVMs'^ and 
in particular through rank-1 POVMs, i.e., projective measure- 
ments'^. Often it is sufficient to consider von Neumann mea- 
surements, but in some in some cases considering more gen- 
eral projective measurements allows for a better minimization, 
an issue still under research (see for instance Ref. "*■*). 
The minimization can be done analytically for some simple 
cases of two-qubits, namely for the class of X states which 
have nonnull entries only on the diagonal and antidiagonal and 
include states with maximally mixed marginals(see Refs4^i^ 
and"*-'"*- for recent developments). On the contrary, the two- 
qutrit case must be handled numerically. 

1. two-qubit states 

Since in phases I and III the density matrix pij corre- 
sponds to an X-state for which vi\f^g^^S{pij\{Bk]) can be 
easily evaluated with a fully analytical way by resorting to the 
method developed in . In this part of the Appendix we give 
a brief review of this method. 

An arbitrary (single-qubit) von Neumann measurement is de- 
fined by a couple of orthogonal projectors _Bo and Bi, which 
can be obtained from |0)(0| and by an arbitrary S'J7(2) 

rotation V: 

Bq = V\Q){Q\V^ Bi = V\l){l\V^ (Al) 

Since V = tl + iy ■ a with + yl + yl + vi = 1. von 
Neumann measurements are parametrized by three indepen- 
dent numbers. 

The key result of is that the minimum of S{pij\{Bk}) 
is always attained for some special values of the parameters 

{tyi+y2y3)^,n= {ty2-yiyi.){tyi-y2yi),k ^ +yl, 
namely 

{fc 0, TO = 0, ri = 0}and 

{fc= 1/2, TO = 0,1/4, n = 0, ±1/8} (A2) 

Therefore the minimization procedure reduces to comparing 
the expressions S{pij\{Bk]) obtained in correspondence of 
these two sets of values. Furthermore, when the two-site re- 
duced density matrix element (p^ )i,4 = 0, which is our case, 
TO and n become irrelevant and S{pij\{Bk]) depends only 
on k. Therefore, we only have to compare S{pij\{Bk]) for 
k = l/2andA: = 0. 

By the formulas in Ref. , for k = 1/2 we have 

Sip,j\{Bk}) = Slip,,) = (A3) 
2 ^062 2 2 '"§2 2 
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where 



^1 = V [(Py)ll - (P*j)33 + (P»j)22 - (Aj)44]2 + 4|(py)23p 

(A4) 

while for fc = we have 

Sip,j\{Bk}) = S2{p^J)= (A5) 
-(1 - Po)i^ log2 - (1 - Po)i^ l0g2 

-p,l^ log, l^-po^^ log, '-^ 



where po = (Pij)ii + (PijOsa and 

^ |(P»j)22 - (P»j-)44| ^ - (P»j)33| ^^g^ 

|(pji)22 + {Pi3)m\ ' + (p^jOasI 

All we must do is take the minimum between (IA3l l and (IA5b : 
inf{B,}5(p„ |{Bfc}) = min{5i, ^2} (A7) 



2. two-qutrit states 

As for the two-qutrit case, we have that the possible von 
Neumann measurements correspond to unitary rotations. 

Bo = V\0){0\V^ ,Bi - V\1){1\VUB2 - V\2){2\V^ 

(A8) 

where now V G SU{3). 

Unfortunately, to proceed forward in the computation of the 



discord, one cannot simply mimic the procedure described 
for qubits. The main difficulty is that no easy, explicit 
parametrization of V G SU{3) by 8 real parameters (the 
group dimension) can be found. We therefore must com- 
pute the discord numerically. Our strategy is to minimize 
S{pij\{Bk}) over a (large) set of randomly-generated uni- 
tary matrices^*. More precisely, we generate a large ensam- 
ble of unitary matrices taken from the uniform distribution 
over the SU{3) group manifold, evaluating S{pij\{Bk}) for 
each matrix. We then keep the minimum as our esteem of 
inf{Bfc} S{pij\{Bk})- To be rigorous, this esteem is to be re- 
garded as an upper bound: however, since we are generating 
a rather large set of random matrices we are confident that the 
bound is very stringent. 

Alternatively, we can use the SU{3) parametrization given 
in Ref.^^. This allows to parameterize SU{3) in terms of 
trigonometric functions of 8 independent parameters, 3 an- 
gles ?7i, ?72, ?73 and 5 phases a, /3, 7, p, a. This parametrization 
makes it apparent that the phases p and cr are completely ir- 
relevant for the computation of the discord, since orthonormal 
projectors (von Neumann measurements) are independent of 
the choice of such phases. This method has the advantage that 
it is based on a more transparent parametrization of von Neu- 
mann measurements. Again, we generate a large ensemble of 
unitary matrices find the minimum of |{i3fe}) 
In all cases under study, the two methods applied led to the 
same results, which provides us with full confidence on their 
reliability. 
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